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ABSTRACT 

We present the AEGIS-X survey, a series of deep Chandra ACIS-I observations of the Extended 
Groth Strip. The survey comprises pointings at 8 separate positions, each with nominal exposure 
200ks, covering a total area of approximately 0.67 deg^ in a strip of length 2 degrees. We describe 
in detail an updated version of our data reduction and point source detection algorithms used to 
analyze these data. A total of 1325 band- merged sources have been found to a Poisson probability 
limit of 4 X 10^^, with hmiting fluxes of 5.3 x 10"^'' erg cm^ s~^ in the soft (0.5-2 keV) band 
and 3.8 x 10^^^ erg cm~^ s~^ in the hard (2-10 keV) band. We present simulations verifying the 
validity of our source detection procedure and showing a very small, < 1.5%, contamination rate 
from spurious sources. Optical/NIR counterparts have been identified from the DEEP2, CFHTLS, 
and Spitzer/IRAC surveys of the same region. Using a likelihood ratio method, we find optical 
counterparts for 76% of our sources, complete to i?^B=24.1, and, of the 66% of the sources that have 
IRAC coverage, 94% have a counterpart to a hmit of 0.9 /xJy at 3.6 jum (m^s=23.8). After accounting 
for (small) positional offsets in the 8 Chandra fields, the astrometric accuracy of the Chandra positions 
is found to be 0'.'8 RMS, however this number depends both on the off- axis angle and the number of 
detected counts for a given source. All the data products described in this paper are made available 
via a public website. 

Subject headings: galaxies: active — galaxies: nuclei — surveys — X-rays: galaxies 



1. INTRODUCTION 

There is intense current interest in the formation of 
galaxies, the history of star formation and accretion 
power in the universe, and the inter-relations between 
these phenomena. This has motivated the investment of 
major observational resources in obtaining images and 
spectroscopy in various areas of the sky. These multi- 
wavelength surveys cover a wide range of areas, depths, 
angular resolutions and wavebands, ranging from very 
deep "pencil" beam surveys in small areas of the sky to 
wider surveys at shallower depths. 

Surveys at X-ray wavelengths provide an important 
component of the multiwavelength arsenal, primarily be- 
cause of the efficiency and sensitivity of X-ray emis- 
sion in selecting Active Galactic Nuclei (AGN). There 
is apparently a strong relati onship between AGN and 
their host galaxy bul ges (e.g. iFerrarese fc Merritti [20001 : 
iGebhardt et al.l [20001 ). so studying the growth of super- 
massive black holes in the context of the evolution of 
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their host galaxies is likely to be a fruitful area of ex- 
ploration. X-ray surveys also enable an estimate of the 
history of accretion power in the universe and the origin 
of the X-ray background (e.g. [Brandt fc HasingeH r2005). 

The deepest X-ray surveys in existence are the Chan- 
dra Deep Fields North and South (hereafter the CDF-N 
and CDF-S, respectively), which each cover an area of 
~ 0. 1 deg^ to nominal depths of ^ 2 — 3 x 10 ~^^ erg cm~^ 
s-i ([Alexander et al.[[2003l: TLuo et al.[[2008[ ). The widest 
surveys to resolve a significant fraction of the X-ray back- 
ground cover of order 10 deg^, but to much sh allower 
(by a factor ^ 100) dept hs (e.g. XBo otes: Murra y et al.l 
[20051 [Renter et al.l [2005': XMM-LSS: [Pim^e et al.l l2004l. 
The major advantage of ultradeep surveys is that they 
are able to probe into the distant universe, to detect 
"typical" objects at high redshift. On the other hand, 
larger areas are required to sample significant large scale 
structures in the universe and hence determine the re- 
lationship between galaxy evolution and local environ- 
mental density. Moreover, larger area surveys are able 
to sample and find unusual, rare objects. For a com- 
plete picture it is clearly also necessary to explore the 
parameter space in between the ultradeep and ultrawide 
surveys. 

Motivated by this, we have obtained deep X-ray obser- 
vations using the Chandra X-ray observatory in a region 
of sky of area ~ 0.5 deg^ known as the Extended Groth 
Strip (EGS), covering the energy range 0.5-7 keV. Ex- 
tensive multiwavelength data in this regio n have been 
obtained as part of the "AEGIS" project ([Davis et al.l 
[2007( ) making it one of the premier datasets to study the 
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co-evolution of black hole accretion and galaxy forma- 
tion. The purpose of the present paper is to describe 
the X-ray dataset and reduction, and present a catalog 
of point sources derived from the Chandra X-ray survey, 
which we designate AEGIS-X. 

The structure of the paper is as follows. In S|2]we de- 
scribe the data and the reduction, including the method 
used to calculate the Chandra point spread function. De- 
tailed descriptions of our source detection and photom- 
etry procedures are given in In 21 the results of the 
source detection and analysis are presented. In Sj5] the 
optical and infrared counterparts to the X-ray sources 
are provided and in ^the conclusions are given. 

2. OBSERVATIONS AND DATA REDUCTION 

The data used in this paper were obtained during two 
Chandra observation cycles. An initial Chandra obser- 
vation in the EGS region was taken in August 2002 
during Cycle 3. These observations of the "Groth- 
Westphal Str i p" (G WS) region have been reported by 
iNandra et all (|2005l . hereafter N05). The majority of the 
AECIS-X data were obtained over the period March- 
December 2005 as part of Chandra Cycle 6. We ana- 
lyze all the AEGIS-X data here in a uniform fashion. 
All the observations were taken using the ACIS-I instru- 
ment without any grating in place. The S2 and S3 chips 
of the ACIS-S array were also sometimes operating dur- 
ing the observations but as these are far off-axis we do 
not consider the data further. Details of the Chandra 
observations are shown in Table [TJ 

2.1. Data reduction 

The data reduction was performed using the CIAO 
data analysis software version 3.3. Basic data reduction 
proceeded in a manner similar to that described in N05. 

Initially, each observation ID (obsID) was analyzed 
separately. We first corrected the raw (level 1) event 
files for any known aspect offsets. Hot pixels and 
cosmic ray afterglows were removed using the CIAO 
acisjrunJiotpix task for EGS observations 1-7, which 
were taken in vfaint mode. Application of this task 
to the EGS-8 data (i.e. the original GWS data), which 
were taken in faint mode, resulted in images with a large 
number of apparently spurious sources, which we identi- 
fied post hoc due to the anomalously low optical/IRAC 
identification rate in this field. Application of the older 
afterglow removal tool acis_detect_af terglow yielded 
source numbers and identification rate similar to the 
other fields. This suggests that the newer CIAO algo- 
rithm fails to identify some afterglows and/or hot pix- 
els in the faint mode data. We therefore adopted the 
acis_detect_af terglow algorithm for the EGS-8 field. 

To produce a level 2 (screened) event file we applied the 
charge-transfer inefficiency (CTI) and time-dependent 
gain correction, removed the ACIS pixel randomization 
and applied PHA randomization, as recommended. In 
the case of vfaint mode data (EGS- 1-7) we also applied 
the ACIS particle background cleaning algorithm. 

To identify periods of anomalously high background, 
which would hamper efficient source detection, we cre- 
ated light curves in the 0.5-7 keV band (excluding 
sources detected using the celldetect algorithm) with 
a bin size of 50s. Backgro und flares were ident ified using 
the procedure described in lNandra et al.l ()2007f ). whereby 



a quiescent background level is determined by calculating 
the count rate limit at which the excess variance of the 
background light curve is equal to zero. For our Chan- 
dra ACIS data, we excluded times when the background 
count rate exceeded 1.4 times this limit (c.f. Nandra 
et al. 2007, who used two times the limit for XMM- 
Newton EPIC-pn background cleaning). This procedure 
was found to be adequate in most cases but some resid- 
ual background flares not identifled by this method were 
removed manually. One observation, obsID 4365 in field 
EGS-8, exhibited a ~25ks period of elevated, but rela- 
tively stable and well-behaved, background which would 
be excluded using the filtering criteria described above. 
As described in N05, the inclusion of this '^25ks increases 
the sensitivity of the observation to point sources and 
therefore it has not been filtered from the data set for 
this work. 

From the cleaned level 2 event file we created images 
in four energy bands, which we designate full band (FB; 
0.5-7 keV), soft band (SB; 0.5-2 keV), hard band (HB: 2- 
7 keV), and ultrahard band (UB; 4-7 keV). We used the 
standard CIAO procedure (merge_all) to produce expo- 
sure maps for each obsID. These account fully for the 
telescope and instrument efficiencies, as well as the chip 
gaps and dithering during the observation. The efficien- 
cies are also a function of photon energy, so creating the 
exposure maps requires an assumption about the energy 
distribution of the detected photons, i.e. the source spec- 
trum. N05 calculated these exposure maps at a single en- 
ergy representative of the detected photons in each band. 
Following this procedure through to the source photom- 
etry stage indicated that these exposure maps produced 
by this method resulted in hard band fluxes that were 
systematically low. We verified this using MARX simu- 
lations (see ^^i^, inputting synthetic sources of known 
flux and running through the entire source detection and 
photometry procedure. Small systematic flux underesti- 
mations were also found for the other bands, but the HB 
suffered the most severe problem, at the ~ 20% level. In 
part, this may be due to an inconsistency between the 
single energy exposure maps, and the F = 1.4 spectrum 
used to convert counts to flux (see i)3.3p . Switching to 
a weighted exposure map based on a F = 1.4 spectrum 
solved the flux problem, and we adopted this for the cal- 
culation of all exposure maps. 

As already discussed, and listed in Table [1] each field 
typically comprises several obsIDs and, furthermore, the 
obsIDs from any given field can overlap the others. It is 
impractical to perform the analysis for the whole EGS 
field, as the resulting images would be too large when 
made at native 0.5 arcsec resolution. We therefore de- 
fined a total of 8 regions (Figure [T|), one for each field, 
which we used to create images and perform source detec- 
tion, eventually merging the source catalogs from these 
regions. 

The field regions were defined based on the limits in 
sky co-ordinates (X, Y) of the first obsID with a bound- 
ary of 20 pixels in all directions. For example, for EGS-1 
the field region is defined by the maximum and mini- 
mum sky co-ordinates for obsID 5841, plus or minus 20 
pixels in both X and Y directions. We then identified 
all other obsIDs which overlap with this field box, re- 
gardless of whether they are nominally part of that field 
(e.g. we would combine overlapping data from EGS- 
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TABLE 1 
Observation Log 



Field^ 


ObsID'^ 




DEC<= 


Start Time<* 


Modc° On timef 


Exposure^ 


Name 




(J2000) 


(J2000) 


(UT) 


(ks) 


(ks) 



EGS- 




5841 


14 


22 


41, 


.88 


+53 


25 


53 


,38 


2005- 


-03- 


■14 


00: 


:04: 


;09 


VFATNT 


44.46 


44.45 


EGS- 


.1 


5842 


14 


22 


41, 


52 


+53 


25 


53 


,09 


2005- 


-03- 


■16 


15: 


:54: 


34 


VFATNT 


46.43 


46.42 


EGS- 


.1 


6210 


14 


22 


43, 


32 


+53 


25 


24 


,26 


2005- 


-10- 


■03 


14: 


:56: 


50 


VFATNT 


46.35 


45.94 


EGS- 


,1 


6211 


14 


22 


43, 


.68 


+53 


25 


26 


,18 


2005- 


-10- 


■12 


11: 


:43: 


;28 


VFATNT 


36.07 


35.64 


EGS- 


,1 


7180 


14 


22 


43, 


.68 


+53 


25 


26 


,18 


2005- 


-10- 


■13 


05: 


:16: 


04 


VFATNT 


20.85 


16.67 


EGS- 


.2 


5843 


14 


21 


32, 


.04 


+53 


13 


42 


,86 


2005- 


-03- 


■19 


17: 


:13: 


09 


VFATNT 


44.47 


44.46 


EGS- 


.2 


5844 


14 


21 


31, 


.68 


+53 


13 


42 


,52 


2005- 


-03- 


■21 


22: 


:37: 


40 


VFATNT 


45.85 


45.85 


EGS- 


.2 


6212 


14 


21 


33, 


.84 


+53 


13 


14, 


.74 


2005- 


-10- 


■04 


22: 


:56: 


06 


VFATNT 


46.29 


46.28 


EGS- 


.2 


6213 


14 


21 


33, 


.84 


+53 


13 


15 


.02 


2005- 


-10- 


■06 


06: 


■ 52; 


;13 


VFATNT 


47.93 


47.46 


EGS- 


.3 


5845 


14 


20 


27, 


.24 


+53 


02 


15 


.70 


2005- 


-03- 


■24 


14: 


:33: 


:31 


VFATNT 


48.41 


48.40 


EGS- 


.3 


5846 


14 


20 


26, 


.88 


+53 


02 


15 


20 


2005- 


-03- 


■27 


04: 


:51: 


;15 


VFATNT 


49.41 


49.40 


EGS- 


.3 


6214 


14 


20 


29, 


.04 


+53 


01 


47 


.16 


2005- 


-09- 


■28 


08: 


:09: 


03 


VFATNT 


47.93 


47.50 


EGS- 


.3 


6215 


14 


20 


29, 


.04 


+53 


01 


47 


.35 


2005- 


-09- 


■29 


15: 


:58: 


09 


VFATNT 


50.30 


48.63 


EGS- 


.4 


5847 


14 


19 


22, 


.44 


+52 


50 


44, 


.27 


2005- 


-04- 


■06 


20: 


:01: 


09 


VFATNT 


44.99 


44.52 


EGS- 


.4 


5848 


14 


19 


22, 


.44 


+52 


50 


44, 


.26 


2005- 


-04- 


■07 


21: 


03; 


59 


VFATNT 


44.46 


44.40 


EGS- 


.4 


6216 


14 


19 


24, 


.60 


+52 


50 


17, 


.38 


2005- 


-09- 


■20 


09: 


35; 


:13 


VFATNT 


49.92 


49.48 




.4 




14 


1 Q 


24 


fin 




ou 


1 7 


.oy 


zuuu- 


flQ 
-uy- 


9*^ 
■zo 




.01. 


.oy 


VFATNT 


AQ Q9 


AQ 

^y.^o 


EGS- 


.5 


5849 


14 


18 


21, 


.60 


+52 


38 


50 


.53 


2005- 


-10- 


■11 


12: 


:47: 


43 


VFATNT 


49.92 


49.46 


EGS- 


.5 


5850 


14 


18 


21, 


.60 


+52 


38 


51 


.16 


2005- 


-10- 


■14 


05: 


:15: 


37 


VFATNT 


45.95 


40.42 


EGS- 


.5 


6218 


14 


18 


21, 


.60 


+52 


38 


49 


.54 


2005- 


-10- 


■07 


05: 


:31: 


;36 


VFATNT 


40.90 


40.58 


EGS- 


.5 


6219 


14 


18 


21, 


.24 


+52 


38 


47, 


.33 


2005- 


-09- 


■25 


15: 


:57: 


04 


VFATNT 


49.90 


49.43 


EGS- 


-6 


5851 


14 


16 


26, 


.04 


+52 


19 


52 


.37 


2005- 


-10- 


■15 


03: 


:03: 


;18 


VFATNT 


36.08 


35.68 


EGS- 


-6 


5852 


14 


16 


26, 


.40 


+52 


20 


06 


21 


2005- 


-12- 


■03 


13: 


:00: 


33 


VFATNT 


10.62 


10.20 


EGS- 


-6 


6220 


14 


16 


24, 


.96 


+52 


19 


46 


.68 


2005- 


-09- 


■13 


09: 


:17: 


01 


VFATNT 


37.90 


34.84 


EGS- 


-6 


6221 


14 


16 


24, 


.96 


+52 


19 


46 


.81 


2005- 


-09- 


■15 


22: 


:11: 


;12 


VFATNT 


4.15 


3.71 


EGS- 


-6 


6391 


14 


16 


24, 


.96 


+52 


19 


46 


.42 


2005- 


-09- 


■16 


20: 


:43: 


01 


VFATNT 


8.72 


7.47 


EGS- 


-6 


7169 


14 


16 


26, 


.40 


+52 


20 


06 


.33 


2005- 


-12- 


■06 


02: 


29; 


;46 


VFATNT 


16.30 


16.03 


EGS- 


■6 


7181 


14 


16 


26, 


.04 


+52 


19 


52, 


.37 


2005- 


■10- 


■15 


21: 


:17: 


;21 


VFATNT 


16.38 


15.98 


EGS- 


■6 


7188 


14 


16 


26, 


.04 


+52 


20 


05, 


.96 


2005- 


■12- 


■05 


04: 


:50: 


:40 


VFATNT 


3.32 


2.58 


EGS- 


■6 


7236 


14 


16 


26, 


.40 


+52 


20 


06, 


.16 


2005- 


■11- 


■30 


19: 


:29: 


:34 


VFATNT 


20.46 


20.02 


EGS- 


■6 


7237 


14 


16 


26, 


.40 


+52 


20 


06, 


.14 


2005- 


■12- 


■04 


05: 


:26: 


;20 


VFATNT 


17.09 


16.93 


EGS- 


■6 


7238 


14 


16 


26, 


.40 


+52 


20 


06, 


.07 


2005- 


■12- 


■03 


10: 


:02: 


:10 


VFATNT 


9.70 


9.53 


EGS- 


■6 


7239 


14 


16 


26, 


.40 


+52 


20 


07, 


.61 


2005- 


■12- 


■11 


08: 


:31: 


:06 


VFATNT 


16.23 


16.03 


EGS- 


■7 


5853 


14 


15 


23, 


.76 


+52 


08 


16, 


.36 


2005- 


■10- 


■16 


20: 


:16: 


;24 


VFATNT 


42.99 


42.47 


EGS- 


■7 


5854 


14 


15 


23, 


.40 


+52 


08 


12, 


.92 


2005- 


■09- 


■30 


23: 


:52: 


;23 


VFATNT 


50.51 


50.07 


EGS- 


■7 


6222 


14 


15 


22, 


.68 


+52 


08 


09, 


.86 


2005- 


■08- 


■28 


17: 


:20: 


;24 


VFATNT 


35.00 


34.69 


EGS- 


■7 


6223 


14 


15 


22, 


.68 


+52 


08 


09, 


.90 


2005- 


■08- 


■31 


05: 


:06: 


:47 


VFATNT 


50.69 


49.51 


EGS- 


■7 


6366 


14 


15 


22, 


.68 


+52 


08 


09, 


.79 


2005- 


■09- 


■03 


06: 


:30: 


;11 


VFATNT 


14.59 


13.88 


EGS- 


■7 


7187 


14 


15 


23, 


.76 


+52 


08 


16, 


.26 


2005- 


■10- 


■17 


19: 


:07: 


:08 


VFATNT 


6.70 


6.59 


EGS- 


■8 


3305 


14 


17 


43, 


.08 


+52 


28 


25, 


.21 


2002- 


■08- 


■11 


21: 


:43: 


;57 


FAINT 


29.41 


29.36 


EGS- 


■8 


4357 


14 


17 


43, 


.08 


+52 


28 


25, 


.21 


2002- 


■08- 


■12 


22: 


:32: 


;00 


FAINT 


86.12 


84.36 


EGS- 


■8 


4365 


14 


17 


43, 


.08 


+52 


28 


25, 


.21 


2002- 


■08- 


■21 


10: 


:56: 


;53 


FAINT 


84.22 


83.75 



^ Field Name: note that EGS-8 field is the original 
ID'^ Nominal position of pointing (J2000)'^ Start date 
time after data screening as described in i|2.1l 



"Groth-Westphal Strip" described by NOS*^ 
and time (UT)^ Observation mode^ Raw exp 



Chandra Observation 
osure time^ Exposure 



2 with EGS-1 where it exists). The exception to this 
is that we did not combine overlapping data with the 
three observations that comprise EGS-8 with the adjoin- 
ing fields (EGS-5 and EGS-6) due to the unfavorable 
combination of small and large off axis angle data and 
therefore data with significantly different point-spread 
function sizes. Source detection tests performed when 
combining data in this way resulted in fewer detections 
compared to considering the EGS-8 data by itself. For 
each overlapping obsID identified we registered the co- 
ordinates relative to the first observation in the stack us- 
ing the acis_align_eveiits task. This performs a source 
detection using the CIAO celldetect algorithm, at a 
3(7 limit, and uses the centroids to the relatively bright 
sources so identified to realign the sky co-ordinates of 
the images. We only performed this realignment in cases 
where the algorithm detects at least four such sources. 
Following the realignment we created a merged event file 
covering the defined field region, and images in the soft, 
full, hard, and ultrahard bands. 
An exposure map for the field region was created by 



summing the exposure maps of the individual obsIDs 
contributing to the stack, which had been created previ- 
ously. A plot of the effective exposure versus survey area 
is shown in Figure [21 comparing AEGIS-X to the CDFs. 
This shows the large increase in area for moderately deep 
exposures (100-200ks) afforded by our survey. 

The final products of this basic data reduction are 
therefore the 8 event files corresponding to each field 
EGS-1 to EGS-8, and the corresponding images and ex- 
posure maps in the four analysis bands. These images 
are of manageable size to perform source detection. They 
also overlap, which help avoid edge effects in the source 
detection, but means we must subsequently delete dupli- 
cate sources detected in more than one nominal field. 

2.2. Point Spread Functions 

To determine accurately the significance and flux of 
our sources requires a good estimate of the Point Spread 
Function (PSF). In the case of Chandra ACIS data this 
depends on the position on the detector and the energy. 
Finding previous methods for estimating the Chandra 
PSF for use in surveys to be unsatisfactory, either be- 
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EG5-1 




24:00.0 22:00.0 14:20:00.0 18:00.0 16:00.0 

Right Ascension {2000) 



Fig. 1. — Mosaic full band image of the Extended Groth Strip showing the location and overlap of the 8 analysis fields. The EGS-8 field 
is the original "Groth- Westphal Strip" data described by N05. 



cause they are not sufficiently accurate (in the case of the 
CIAO task mkpsf ) or because they are prohibitively slow 
(in the case of ChaRT, the Chandra ray-tracer simulator) 
we decided to develop our own procedure for calculating 
the PSF. The most important difference in the analysis 
of the data in this paper compared to that described in 
N05 is the use of the new calculations of the PSF. 

N05 used the CIAO task mkpsf to estimate the PSF, 
which calculates a PSF image based on a look-up table 
based on calculations u sing the Chandra Ray Tracer PSF 
simulator ChaRT ( Carter et al .l [2003.1 . ChaRT simulates 
photons passing through the best available model of the 
High Resolution Mirror As sembly (HRMA) , b ased on ex- 
tensive pre-launch testing (jJerius et al.lll995h . which in- 
cludes details of the support structures, stray light baf- 
fles, and an independent model of each of the mirror 
segments, and thus provides the best modelling of the 



PSF for any energy or position in the focal (detector) 
plane. The available ChaRT-generated library files for 
mkpsf cover a relatively small number of positions and 
energies, and the tool must therefore interpolate to gen- 
erate the necessary PSFs. In addition, for large off-axis 
angles the PSF images created by mkpsf have the edges 
clipped, as the PSF can become larger than the gener- 
ated image. 

We have therefore taken an alternative approach, 
and ca lculated the PSF using the MARX simulator 
(jWise et a l. 20031. MARX is the detector simulator, but 
is also able to simulate photons passing through a simpli- 
fied version of the ChaRT HRMA model. In particular, 
the MARX HRMA model does not include the physical 
support structure and thus PSFs generated by MARX 
lack some of the complicated substructure seen in the 
ChaRT PSFs. However, knowledge of the detailed struc- 



AEGIS-X 



5 



2500 




r , , ij 

10 100 1000 

Effective exposure (ks) 

Fig. 2. — Effective exposure time versus survey area for AEGIS- 
X, the Chandra Deep Fields North and South and the Extended 
Chandra Deep Field South. These were calculated by taking the 
value of the full band exposure map at each pixel and dividing by 
the on-axis effective area. The CDFs reach much deeper exposure 
but AECIS-X provides a large increase in area compared to these 
surveys at moderate depth. 

ture within the PSF is not needed for our point source 
detection and photometry methods. We only require an 
aperture which contains a particular Enclosed Energy 
Fraction (EEF), and thus the simplified mirror model is 
sufficient for our needs, as well as requiring less computer 
processing time. 

We have generated look-up tables of the PSF for a 
range of EEFs over a fine, evenly spaced grid (10 pixel 
spacing) covering the ACIS-I detector. At each position, 
an image of the PSF was generated with MARX for a 
monochromatic source initially of energy IkeV (repre- 
sentative of our soft band), using 200,000 input rays. 
The effects of quantum efficiency and the filter transmis- 
sions were turned off, as these properties do not affect the 
PSF, thus maximizing the number of detected photons 
for a given number of input rays, improving the count- 
ing statistics in the PSF image without requiring addi- 
tional ray traces and longer processing time. At each po- 
sition, counts were extracted within elliptical apertures 
(found to best describe the overall shape of the PSF) 
with increasing semi-major axes. For each value of the 
semi-major axis, we calculated the semi-minor axis and 
orientation angle of the ellipse from the moments of all 
counts within a circle of radius equal to the semi-major 
axis. We calculated the EEF for each of the set of ellip- 
tical apertures, and used interpolation to calculate the 
ellipse parameters for an aperture containing a particu- 
lar EEF (specifically 50, 60, 70, 80, 90, and 95%). The 
total counts (for normalization of the EEF) were taken 
from the entire area of the detector. When sources fell 
close to the edge of the detector chips or the chip gaps, 
MARX was altered to shift the nominal position of the 
detector in the focal plane by around 100 pixels so the 
PSF was fully sampled. The ellipse parameters for the 
desired EEFs were stored in look-up tables, providing a 
very fast method of determining the PSF, using the clos- 
est grid position. The PSF varies little on our sampling 
scale of 10 pixels, so we do not need to interpolate. The 



procedure was repeated using representative monochro- 
matic energies for each of our energy bands (FB: 2.5 keV, 
HB: 4 keV, UB: 5.5keV). 

The PSF is defined at a position on the detector (fixed 
relative to the mirror). However, our X-ray data are the 
result of a number of observations (obsIDs) with differ- 
ent orientations and pointing directions, merged to create 
images with the maximum possible exposure. The source 
detection and photometry described below requires ex- 
traction of counts using a single, circular aperture from 
stacked images which have different PSFs. We therefore 
calculated exposure weighted PSFs for each candidate 
source in the merged images. The position in each of 
the constituent obsIDs was found, and in each obsID the 
average value of the exposure map at that position (ex- 
tracted over approximately 80 surrounding pixels) was 
calculated, and the PSF was retrieved from the MARX 
table. Each of the MARX PSFs were converted from el- 
lipses to circles with a radius equal to the square root of 
the product of the major and minor ellipse axes. For each 
candidate source the PSFs were then combined, weight- 
ing each by the ratio of exposure in a particular obsID 
to the exposure in the merged image. This was repeated 
for each energy band. 

3. SOURCE DETECTION AND PHOTOMETRY 

The point-source catalog was created using an updated 
version o f the source detect ion algorithm described by 
N05 and iLaird et all (|2006l ). Briefly, for the AEGIS- 
X data, candidate sources were initially identified in 
each individual field in a number of different energy 
ranges, their significance computed, and a threshold 
applied. The sources considered significant were then 
band-merged, and photometry performed to estimate 
the fluxes in several energy ranges. Finally the source 
catalogs for the individual (overlapping) pointings were 
merged to remove duplicate sources for the final catalog. 
These steps are described in detail below. 

3.1. Source detection 

The detection algorithm described in N05 involves pre- 
dctection at a low significance threshold followed by aper- 
ture extraction of the counts to determine the source 
significance. The pre-detection is performed using the 
CIAO wavelet-detection algorithm wavdetect, run on 
the unbinned, stacked images, at wavelet scales of 1, 
\/2, 2, 2^2, 4, 4V2, 8, 8^2, and 16 pixels. Only po- 
sitions with at least 10% of the maximum field exposure 
were considered, and the procedure is run at a signifi- 
cance threshold of 10"'' on images in all four bands. The 
minimum exposure requirement means that in practise 
our source catalog will not cover the full 0.64 degree^ 
of the survey: the least sensitive 3.2% of the total area 
is excluded. The chosen low threshold for wavdetect is 
likely to result in a large fraction of spurious sources, but 
very few real sources should be missed. To determine the 
source significance we first calculated the 70% EEF PSFs 
for the stacked image, by exposure weighting the individ- 
ual PSFs from each individual exposure in the stack as 
described above. 

Using the 70% EEF PSF (determined a priori to be the 
most efficient radius for source detection) we extracted 
the total counts for each candidate source, and the effec- 
tive exposure, the latter being the value of the exposure 
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map averaged over all the pixels in the extraction region. 
Background counts were determined by first masking the 
image to remove pixels within the 95% EEF of each can- 
didate wavdetect source. We then extracted background 
counts from an annulus with an inner radius equal to the 
1.5 times the 95% EEF radius for the source, and an 
outer radius 100 pixels more than this. An average ex- 
posure map value was calculated for this background area 
also. The counts in the background area were than scaled 
to the source region by the ratio of the source and back- 
ground areas, and the ratio of source and background 
average exposure. 

For each candidate source in each energy band, the 
Poisson probability that the total counts in the source 
region would be observed based on the expected back- 
ground counts was then calculated. Following N05, we 
adopted a significance of 4 x 10~^. Typically, when run 
at a significance of 10~^ wavdetect identifies 400-500 
candidate sources in each AEGIS-X field, but more than 
half of these arc found to be insignificant - and hence 
probably background fluctuations - when we perform the 
aperture extraction and Poisson probability calculations. 
Because we initially masked out all the wavdetect can- 
didate sources we may therefore have incorrectly masked 
out the largest positive background fluctuations, leading 
to an underestimate of the true background level. For 
this reason we perform a second pass of the source de- 
tection, masking out only the sources with Poisson prob- 
ability < 4 X 10^^ when calculating the background. We 
then recompute these probabilities with the new back- 
ground estimates and reapply the threshold. This sec- 
ond calculation of the Poisson proba bility is a refinemen t 
of the method described in N05 and iLaird et al.l (|2006D . 
Technically further iterations of this kind could be ap- 
plied but in practice there is negligible difference in the 
background estimates after two iterations. 

When run at low thresholds, such as we have used for 
the seed catalog, wavdetect sometimes detects the same 
source twice, particularly at larger off-axis angles where 
the counts are spread out over a wide area. Visually 
inspecting the images it is clear that only one source is 
present. During the source detection phase we therefore 
check for sources that are separated by less than 5 pixels. 
If there are any then we remove from the list the source 
position with fewer counts in the 70% EEF area. There 
are a total of 4 such cases in AEGIS-X. 

3.2. Band-merging 

Thus far the source detection procedure has been car- 
ried out separately for the soft, full, hard and ultrahard 
band images. While it is possible that a source will be 
detected in just one of these bands, in practice the ma- 
jority of course are detected in several, necessitating a 
band-merging procedure to produce a single source cata- 
log for that field. In practice most objects are detected in 
the full band, so we matched the soft, hard and ultrahard 
band source lists with the full band. The positional ac- 
curacy of Chandra depends on the off-axis angle (OAA), 
due the position dependence of the PSF, so we adopted a 
variable radius for source matching. The cross-matching 
tolerances are shown in Table [21 and are based on posi- 
tional accuracies as determined in i i4.4l We used a ia 
radius, combining the uncertainty in the two positions 
in quadrature. If more than one cross-band counterpart 



TABLE 2 

Cross-band matching radii 



OAA^ 


' IcT 


Tolerance'^ 


(arcmin) 


(arcsec) 


(arcsec) 


0-3 


0.30 


1.30 


3-6 


0.57 


2.44 


6-9 


1.13 


4.79 


> 9 


1.67 


7.08 



^ Off-axis angle*^ Icr posi- 
tional uncertainty determined 
from MARX simulations (see 
text)"^ Cross-band match toler- 
ance 



was identified within the tolerance the closer match was 
selected. There was only one case of this in the survey, 
in EGS-5. Two soft band sources are detected within the 
match tolerance radius from a single full band source. Vi- 
sual inspection of this apparent "double" suggested that 
it was in fact a single object detected twice with a slightly 
different position. The soft position with the larger off- 
set from the full band source was therefore removed from 
the catalog. 

The source positions used in the band merged catalog 
are decided as follows: if a significant full band counter- 
part exists then the full band position is used. If not 
and a soft band counterpart exists then the soft band 
position is used. The hard band position is used for hard 
band sources without full or soft counterparts. Ultrahard 
band positions are used for sources only detected in this 
band. 

After performing the band merging we visually in- 
spected the images and source list in each of the fields and 
checked that the correct cross-band counterparts were 
identified. 

3.3. Photometry 

As with source detection, photometry was carried out 
on each of the eight individual EGS fields. For each 
source in the band-merged catalog we extracted total 
counts and effective exposure values from the 90% EEF 
area in each analysis band. Background counts and ex- 
posure values were then extracted in these same bands 
from an image in which the 95% EEF area of all the sig- 
nificantly detected sources (regardless of band) had been 
masked out. The background extraction region was de- 
fined in the same way as it was for the detection ( i)3.ip . 
Background counts were once more scaled based on the 
ratio of the areas of the source and background cells and 
the ratio of the average effective exposures in all pixels 
in the cells. These scaled background counts were then 
subtracted from the total counts to give the net source 
count. 

To identify sources that were close to each other, and 
hence possibly have confused photometry, we identified 
those that were separated by more than the match toler- 
ances given in Table[2]but whose 90% EEF radii overlap. 
These are fiagged in the catalog as having photometry 
that is likely contaminated by a nearby object. In some 
cases, the actual position can fall within the 90% EEF 
radius of the nearby source. In these cases the photom- 
etry will be heavily contaminated by the nearby source, 
and they are flagged as being confused. In such cases the 
positions may also be highly inaccurate. 
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We calculated source fluxes and la confidence lim- 
its from the detected counts using a Bayesian method 
that correc t s for the Eddington bias and is based on 
iKraft et al.l()199lD . This methodolo gy assumes the range 
of possible source fluxes is a continuous distribution fully 
described by the observed data and prior knowledge of 
the physical system. We apply Bayes' theorem, 



P{S\N, B) cx C{N\S,B)tt{S) 



(1) 



where P{S\N^ B) is the posterior distribution function 
for the source rate, S. £{N\S,B) is the Poisson like- 
lihood of obtaining the observed data, N counts, for a 
given source rate, S, and background B. Thus, 



C{N\S,B) = 



jS + B ) 
Nl 



N 



(2) 



This approach naturally incorporates the contribution 
of both source and background components to the to- 
tal observed counts, and thus provides a better descrip- 
tion of the nature of faint sourc es than classical ap- 
proaches norm ally emp loyed fe.g. iGiacconi et alJ 1200^ 
[Alexander et a l. 2003; Kim et al.il2007[ ). ttIS) is the prior 
probability distribution, and reflects our prior knowl- 
edge o f the range of possible source fluxes. iKraft etall 
()1991h adopted a constant uninformative prior, requir- 
ing only that the source flux must be greater than 0. 
However, the distribution of fluxes of X-ray sources (the 
logTV— logs') are found to be well described by a bro- 
ken power-law (iGeorgakakis et aTl 20081). This intro- 
duces Eddington bias (jEddingtonl [1940) which affects 
the classical estimate of source fluxes since a larger num- 
ber of faint sources are scattered to higher fluxes than 
vice- versa. Thus the fluxes of faint sources are generally 
over-estimated. We can correct for the Eddington bias 
with our Bayesian method, by adopting a prior based on 
the observed logTV— logS relation. To approximate the 
logiV— logS we use a faint end slope of -1.5 and a bright 
end slope of -2.5 for each of the analysis bands. We use 
break fluxes for each band as given in Georgakakis et aTl 
(|2008l Table 2). The best estimate of the source flux is 
then obtained by finding the mode of the posterior distri- 
bution. Differentiating equation [1] provides an analytic 
solution for a prior with power-law slope (3, 



S 



^(^N-B + f3+ ^{N ~B + I3f+ ABI3) (3) 

This expression is easily extended to the case of a broken 
power-law. An additional condition for the existence of 
a solution is introduced. 



{N - B + I3f> -ABl3 



(4) 



When there is no solution, we only quote upper limits on 
the source fluxes. 

For the Bayesian approach, confidence limits are ob- 
tained by integrating the posterior distribution. 



CL= / C{N\S,B)Tr{S) dS 



(5) 



^" (S + b)^ 



e-(^+''V(S) dS 



where CL is a given confidence level, and if is a normal- 
ization constant. For a faint-end slope of (3i — —1.5 the 



posterior distribution diverges at the faintest fluxes. Ex- 
tending the power-law distribution to the faintest fluxes 
is clearly unphysical, and makes normalization impossi- 
ble. We thus set a lower limit to the possible source fluxes 
which is 10 times fainter than the detection limit at the 
position of each so urce, given the background. Following 
IKraft et al.l (|1991h . we adopt Highest Posterior Density 
intervals, which minimize the size of the confidence inter- 
val for the given confidence level (68%), and thus sets the 
tightest limits on the source flux. These are calculated 
by numerical integration of equation 5. 

We use this method to calculate the fluxes of all sources 
in each of the four bands, regardless of their Poisson 
probability or net counts. While many sources are not 
detected in all four analysis bands, we can nonetheless 
often determine reasonably reliable fluxes in those un- 
detected bands. Effective on-axis source count rates 
were calculated by dividing the net counts with the 
average value of the exposure map, and aperture cor- 
recting. Count rates were converted to fluxes using a 
r = 1.4 spectrum with G alactic iVn of 1.3 x 10^" cm~^ 
(jDickev fc Lockmanlll990l ). The count rates in the full, 
hard and ultrahard bands were also extrapolated to the 
standard bands: 0.5-10, 2-10, and 5-10 keV, respec- 
tively. 

For comparison, we also calculated fluxes using a clas- 
sical method, with the la confidence ranges on the de- 
te cted counts de termined according to the prescription 
of lGehrelsl ([l98l . In Figure [3] we compare the Bayesian 
and classical methods for calculating fluxes, for full band 
sources. For bright sources, /0.5-10 kcV ^ lO^^'' erg 
cm~^ s~^, the fluxes and errors from both methods agree 
well. However, the fluxes of the faintest sources do not 
agree because the classical method fails to correct for the 
Eddington bias or accurately describe the errors. In such 
cases the classical method can over predict the fluxes by 
up to factor of 2. 

3.4. Hardness ratios 

HRs were calacu lated using a Bayesian approach, fol- 
lowing IParF^eF^O ([2OO6;) . Their method models the de- 
tected counts as a Poisson distribution and gives error 
bars and reliable HRs for sources with both low and high 
counts. To calculate the HRs with the Bayesi an approach 
we u se the BEHR^ package (ver. 11-08-2007. iPark et all 
I2OO6D . For this we use a flat, non-informative prior and 
use the Gaussian quadrature method for sources with less 
than 20 net counts in either the soft or hard bands, and 
the Gibbs sampler method for all other sources. Effective 
area corrections to on-axis values are input to BEHR to al- 
low on-axis HRs to be calculated and we use the mean 
of the posterior p robability distrib ution as the best esti- 
mate for the HR (jPark et al.ll2006D . 

For comparison, we also calculated HRs using the clas- 
sical method: HR= {H - S)/{H + S), where H and S 
are simply hard and soft band net counts respectively, 
corrected to on-axis values. Where a source has only 
upper limits in the soft or hard bands, the HR is desig- 
nated as -1-1 and —1, respectively. Figure [4] compares the 
two methods of calculating HR. There is good agreement 
for the bright, highly significant sources, but for weak 
sources the two methods disagree, particularly at low 

See |http: / /hea-www. harvard.edu/ AstroStat/BEHR/l 
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Fig. 3. — Comparison of a classical and a Bayesian method 
(which corrects for the Eddington bias) of calculating fluxes. (Top) 
Classical full band flux (0.5—10 keV) versus Bayesian flux. Errors 
bars show the la confldence limits for both methods. (Bottom) 
Classical full band flux versus the ratio of Bayesian to classical 
flux. At the brightest fluxes both methods show good agreement. 
At fainter fluxes the classical method can over predict the fluxes 
by up to 100%, due to the Eddington bias. The same trend is seen 
in the soft, hard, and ultrahard bands. 

and h igh H R values, c o nfirm ing the results of lPark et al.l 
(|200^ and lKim et all (|200l . 

3.5. Sensitivity maps 

For a number of applications, one needs to know the 
flux sensitivity of the images as a function of position and 
energy band. We have calculated sensitivity maps, and 
hence also the area of the AEGIS-X s urvey as a function 
of flux, using the method described bv lGeorgakakis et al.l 
(|2008). This accurately estimates the probability of de- 
tecting a source with a given X-ray flux accounting for 
observational effects, flux estimation biases, and the frac- 
tion of spurious sources expected in any source catalog. 
The backbone of the method is the Poisson probability 
distribution. For a detection cell with local background, 



1.0 F 




-1.0 -0.5 0.0 0.5 1.0 
Classical hardness ratio 

Fig. 4. — Comparison of a classical with a Bayesian method 
of calculating HR for all sources in the AEGIS-X catalog. The 
red points are sources with Poisson probability p < 10~* in both 
soft and hard bands, for which the classical and Bayesian methods 
agree. 

B, one can estimate the minimum integer number of total 
counts, L, for which the Poisson probability that fluctua- 
tions of the background produce at least L counts is less 
than the adopted threshold for source detection, 4x 10~^. 
Repeating this exercise for different cells across the im- 
age one can determine L as function of position on the 
detector. The 2-D image of L is the sensitivity map. The 
sensitivity map determined in this way does not require 
any assumptions on the spectral shape of the source. A 
useful 1-D representation of the sensitivity map is the 
total detector area in which a source with flux fx can be 
detected. The sensitivity curve is estimated by adopt- 
ing a Bayesian approach. For a source with flux fx and 
a given spectral shape (F — 1.4 in this paper) we de- 
termine the probability of detection in a cell with mean 
background B and detection limit L. The total observed 
counts in the cell are T = B + S, where S is the mean ex- 
pected source contribution. In practise this depends on 
the observation exposure time, the vignetting of the field 
at the position of the cell, and the fraction of the total 
source counts in the cell because of the PSF size. The 
probability that T exceeds L is given by the cumulative 
Poisson probability. 

Again, for comparison purposes, we have also calcu- 
lated the sensitivity maps using a "classical" method (as, 
e.g. N05), whereby we assign a single limiting flux to each 
detection cell. The sensiti vity curves in diffe r ent en ergy 
bands estimated using the iGeorgakakis et all (|2008D and 
classical methods are shown in Figure \5\ In the classical 
method there is a limiting flux below which the se nsitiv- 
ity dr ops to zero. In contrast, in the Georgakakis et al.l 
(|2008f l method, at faint fluxes as — > the Poisson prob- 
ability converges to pthresh = 4 x 10~^, because this is 
the finite probability of a random background fluctuation 
above the limit L. 

4. POINT SOURCE CATALOG 
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Fig. 5. — Sensitivity curves in the full, soft, hard, and ultra- 
hard bands. The soli d lines a r e the area curves using the Bayesian 
method of Georgakakis et al.l I I2008I) described in i|3.5l The dashed 
line shows the classical area curves that do not account for the Pois- 
son probability density distribution of the total observed counts in 
the detection cell. 



TABLE 3 
Positional offsets 

APPLIED to individual 
FIELDS 



Field 


ARA 


ADEC 




(arcscc) 


(arcsec) 


EGS-1 


-0.081 


0.046 


EGS-2 


-0.031 


0.060 


EGS-3 


0.025 


0.048 


EGS-4 


-0.151 


0.086 


EGS-5 


-0.197 


-0.004 


EGS-6 


-0.068 


-0.051 


EGS-7 


-0.102 


-0.002 


EGS-8 


-0.072 


0.059 



To produce the final catalog we need both to remove 
duplicate sources detected in more than one of the (over- 
lapping) fields, and to register all the images to the same 
astrometric frame. To do so, we first matched the X- 
ray catalogs to t he DEEP2 optical photometry catalog 
(|Coil et al.ll2004D and the median RA and DEC offsets 
were calculated. The X-ray positions were shifted to be 
in agr eement with the DEEP2 images (i.e. the SDSS 
frame: [Coirerall[200l using the offsets given in Table [3l 

The source lists were then combined by finding the 
closest source within a match radius of 5". Typically the 
offsets were much smaller than this (see HAA^i . however 
as sources are matched only between individual fields and 
not within them there is no danger of close sources being 
erroneously combined. We visually compared the com- 
bined source list to the individual field source lists to con- 



TABLE 4 

Number of sources detected using different methods 



Method 


Pthresh 


FB^ 


SB'' 


HB^ 


UB^' 


Merged 


This work 
wavdetect 


4 X 10-^ 
1 X 10-'^ 


1221 
1160 


1032 
968 


741 
696 


350 
311 


1325 
1258 



FB = full band (0.5-7 keV), SB = soft band (0.5-2 
keV), HB = hard band (2-7 keV) and UB = ultrahard 
band (4-7 keV). 



TABLE 5 

Sources detected in one band but not 

ANOTHER 



Detection 




Non-detection band 


band (keV) 


Full 


Soft Hard Ultrahard 


Full (0.5-7) 




274 497 876 


Soft (0.5-2) 


85 


466 748 


Hard (2-7 ) 


17 


175 ... 396 


Ultrahard (4-7) 


5 


66 5 



firm this. Where a source is detected in more than one 
field, we quote the source properties (e.g counts, proba- 
bilities, fluxes) for the field in which the source has the 
smaller OAA relative to the field center. Usually the 
source properties will be very similar in the two overlap- 
ping frames in which it is detected, but can differ subtly 
as e.g. the backgrounds will in general be different. The 
original detected field can be identified by the 'FieldJD' 
tag in the source catalog. The final catalog was ordered 
and numbered according to increasing RA. 

Once double sources are removed there are a total of 
1325 distinct sources in the final catalog (Table [9|). De- 
tails of how many sources are detected in each individual 
band are given in Tableland a summary of the number 
of sources detected in one band but not another is given 
in Table O A large number of sources are detected (at 
p < 4 X 10~^) in just one band, specifically, 115, 85, 14, 
and 2 in the full, soft, hard, and ultrahard bands respec- 
tively. The brightest significant source in the catalog has 
a total of 4290.5 net counts in the FB and the faintest 
has 4.2 net counts. 

4.1. Flux limit 

Using the sensitivity curves in Figure [5] we determine 
the limiting fluxes in each band, defined as the flux to 
which at least 1% of the survey area is sensitive, to be 
2.37 X 10-16 (FB; 0.5-10 keV), 5.31 x 10"!^ (SB; 0.5- 
2 keV), 3.76 x 10'^^ (HB; 2-10 keV), and 6.24 x lO^i^ 
erg cm"^ s~^ (UB; 5-10 keV). The total area of the sur- 
vey is 0.64 degree^. Using classical sensitivity curves 
which do not correct for the Eddington bias, as in N05, 
gives limiting fluxes that are approximately a factor of 
2 greater. Because the sensitivity varies dramatically 
across the images, we also give the flux limits corre- 
sponding to 50% and 90% completeness in Table [H Fig- 
ure [5] shows the distribution of X-ray flux in each of 
the four bands. The fluxes cover a range of over three 
dex, with ~50% of detected soft and hard band sources 
having fluxes of less than 9.7 x 10^^^ and 4.5 x 10^^^ 
erg cm"^ s~^, respectively. 

4.2. Comparison with wavdetect 
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Fig. 6. — The distribution of fluxes for sources detected in the 
full, soft, hard, and ultrahards bands, using the Bayesian calcula- 
tion of the fluxes. 



TABLE 6 

AEGIS-X FLUX LIMITS 





Completeness limit^ 


Band 


1%^ 


50%'= 90%'= 


Full 


2.37 


13.03 40.27 


Soft 


0.53 


3.35 11.09 


Hard 


3.76 


20.65 62.37 


Ultrahard 


6.24 


38.50 124.45 



Flux to which 1, 50 and 
90% of the survey area is 
complete. '= Fluxes are in units of 
10"^^ erg cm"^ s'V 

The standard procedure for source detection with 
Chandra A CIS data is t he wavelet-based algorithm 
wavdetect (jFreeman et al.l f2002'). In order to compare 
our detection method, and verify that it is correctly iden- 
tifying significant sources, we also performed source de- 
tection using the CIAO wavdetect algorithm in each 
of the 8 fields. Running wavdetect with a detection 
threshold of 10~^, c ommonly used fo r ACIS-I s urveys 



Alexander et al.ll2003l : lYang et al.|[2004 : Wang et all 



|2007| ). and creating a band-merged catalog for the whole 
of the EGS we found the total number of independent 
sources detected was 1258, compared to 1325 using our 
method, and less sources were detected in every band 
(Table HI . 

Matching the two catalogs showed that there were 1208 
sources in common. There were 59 sources that were only 
in the wavdetect catalog and 122 sources that were only 
in the ^i?G'/5'-X catalog. To assess the likelihood of these 
sources being real we searched for optical and infrared 
(IR) counterparts in the DEEP2 and Spitzer/IRAC pho- 



tometry catalogs covering the EGS (see ^ for details 
of the optical and IR data and the method for identify- 
ing secure counterparts). Real X-ray sources are likely 
to have a higher counterpart identification rate. We find 
secure optical counterparts for 59% of the AEGIS-X-on\y 
sources, compared to an optical identification rate of 47% 
for the wavdetect-only sources. Of the 73 AEGIS-X- 
only sources that are covered by the Spitzer/IRAC data 
88% have secure IR counterparts, compared to 40% of the 
43 wavdetect-only sources covered by the data. There- 
fore while some fraction of the 59 X-ray sources that are 
in the wavdetect catalog, but have been missed by our 
source detection, may indeed be real, we can be confident 
that the extra sources included in the AEGIS-X catalog 
are at least more likely to be real X-ray sources than 
those that were missed. 

4.3. Comparison with N05 GWS analysis 

The three obsIDs that constitute EGS-8 were pre- 
viously reduced and analyzed as part of the "Groth- 
Westphal Strip" survey and a point source catalog was 
presented in N05. In this more recent reduction and anal- 
ysis of the data there are a number of differences in the 
source detection results, compared to the N05 results. As 
was stated in ^ 12. II this field was analyzed independently 
of the rest of the EGS data and therefore includes the 
same data as in our previous work. In both analyzes 158 
independent sources were detected in the field, however 
the IDs of all of the sources are not the same. There 
are 10 new AEGIS-X sources that were not detected in 
N05 and 9 sources in the N05 GWS catalog that are not 
included in the ^£'G/S'-X catalog (Tablel?]). One further 
GWS source (cl3) was not detected in EGS-8 but was 
detected in EGS-6, which overlaps with the field, and 
is therefore included in the final source list. The high 
op tical identification r ate o f the sources, as determined 
bv iGeorgakakis et al.l (|2006D . suggests that the majority 
of the 9 sources are probably real. However, in order to 
maintain a well controlled sample we have not included 
them in our final AEGIS-X catalog. 

There are several small differences in our reduction and 
analysis that could lead to the different source lists. The 
EGS-8 data used in this work has a longer total exposure 
(^ 7ks) and the PSFs were calculated using a different 
method and are on average 7% larger than before. How- 
ever the probable main reason that the sources are not 
detected in this analysis is the refinement of our source 
detection method from the method used in N05. The 
second pass of source detection described in ^ 33. II leads to 
25% higher background levels than in N05 and therefore 
faint sources with Poisson probabilities close to the de- 
tection threshold, such as those in Table [71 will no longer 
pass the detection criteria. The Poisson probabilities of 
the 9 N05 sources using the new higher background are 
shown in Table [Tj several sources were just below our de- 
tection criteria for being included in the AEGIS-X cata- 
log. 

4.4. Astrometric accuracy 

The astrometric accuracy of the Chandra positions 
was estimated using MARX to generate simulated fields. 
Synthetic point sources were added to the simulated im- 
ages adopting a single power-law log N — log S distribu- 
tion with parameters tuned to reproduce the faint end of 
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TABLE 7 

GWS SOURCES FROM N05 NOT INCLUDED IN AEGIS-X CATALOG 



N05 


RA J2000'' 


Dec J2000'' 






Det.'' 


Rab" 


Comments^ 


„ . b = d 




cat. no. 


(J2000) 


(J2000) 






bands 


(mag) 








cl4 


14 16 59.26 


+52 34 36.04 


1.3 X 10- 


-6 


fh 


25.03 




5.5 X 10- 


6 


c35 


14 17 18.89 


+52 27 43.74 


1.1 X 10" 


-6 


f 


23.47 


chip gap 


4.1 X 10" 


■6 


c45 


14 17 25.28 


+52 35 12.08 


4.0 X 10" 


-6 


f 


>26.0 




1.8 X 10" 


■4 


c68 


14 17 39.06 


+52 28 43.78 


3.2 X 10" 


-7 


fs 


24.92 


chip gap 


2.3 X 10" 


■3 


c69 


14 17 39.31 


+52 28 50.16 


6.3 X 10" 


■8 


fh 


23.26 


chip gap 


7.6 X 10" 


■4 


c73 


14 17 42.86 


+52 22 35.21 


1.9 X 10" 


■8 


f 


>26.0 


chip gap 


3.1 X 10" 


3 


c80 


14 17 47.01 


+52 25 12.07 


1.2 X 10" 


-6 


£h 


22.54 




1.3 X 10" 


■5 


cl02 


14 17 56.92 


+52 31 18.47 


3.7 X 10" 


-6 


f 


25.10 




3.6 X 10" 


■5 


cl25 


14 18 08.06 


+52 27 50.36 


2.4 X 10" 


6 


f 


>26.0 




2.7 X 10" 


■4 



^ Values from N05.'' Lowest fal se detection probability found for the four analysis 
bands. Optical identification from iGeorgakakis et al.l ||2006D .'^ Values from this work. 



TABLE 8 

Absolute astrometric accuracy of sources 
IN THE AEGIS-X catalog 



OAA 
(arcmin) 


Net full band counts 
5 ^ 50 50 < 5 ^ 100 S > 100 
(arcsec) (arcsec) (arcsec) 


0-4 


0.56 


0.36 


0.14 


4-8 


0.96 


0.61 


0.44 


> 8 


1.33 


0.87 


0.57 



the 2-10 keV X-ray number counts. The slope of the dif- 
ferential counts was set to /? = —1.58 and the normahza- 
tion is fixed so that there are 7000 sources per square de- 
gree brighter than /x(2 — lOkeV) — 10"^^ ergs~^ cm~^, 
i.e. similar to the observed density of X-ray sources 
(|Kim et all [20071 : IGeorgakakis et al.H2008f) . This choice 
of parameters will over predict the number of bright 
sources but for the purposes of determining the positional 
error this is irrelevant. We adopt an X-ray spectrum 
with r = 1.4 and a total exposure time of 200ks. Sim- 
ulated ACIS-I event files are constructed by randomly 
placing within the Chandra field of view point sources 
with fiuxes in the range fx {2 - lOkeV) = 5 x 10"" - 
10~^^ ergs""'^ cm~^. MARX does not simulate the Chan- 
dra background. This is added to the simulated images 
using the quiescent background event files produced by 
the ACIS calibration team using blank sky observations. 
A total of 10 such simulations were performed. Each 
simulated field was ran through the entire source detec- 
tion procedure. In each simulation, there are on average 
about 180 sources with detection probability < 4 x 10~^ 
in at least one of the four analysis bands. The detected 
source positions were then compared to the input posi- 
tions of the sources in the simulation. It should be noted 
that the astrometric accuracy determined this way, us- 
ing simple, single field images, accounts only for statis- 
tical errors and does not take into account any system- 
atic errors that may be present using the more complex 
AEGIS-X data. 

As expected, the positional accuracy depends on both 
the OAA, and the number of source counts. How- 
ever the simulations showed only a mild dependence on 
source counts and therefore for the purpose of cross- 
band matching we chose radii that depended only on 
OAA. The estimated la positional accuracies are shown 
in the second column of Tabled averaged over all source 
counts. 



The absolute astrometric accuracy of the Chandra po- 
sitions in the ^iJG/S'-X catalog were verified by matching 
the AEGIS-X sources to the DEEP2 optical data from 
the AEGIS project (see SjST] for details), for which the 
absolute positional uncertainty is 0'.'2 RMS (A. L. Coil, 
priv. communication). As described above, to combine 
the catalogs from each of the eight fields we first reg- 
istered the positions to match the astrometric frame of 
the DEEP2 catalog using the offsets given in Table [H 
Matching the final catalog to DEEP2, and considering 
only secure optical counterparts, we find that the me- 
dian and RMS positional offsets of 868 matched sources 
are 0'.'40 and 0'.'79, respectively. The dependence of the 
astrometric accuracy on off-axis angle and source counts 
are illustrated in Figure [71 We determine the la astro- 
metric accuracy of sources in the AEGIS-X catalog by 
breaking the sample down into nine bins based on OAA 
and net counts (Table [8|) . 

4.5. False source contamination 

The number of false sources in our catalog was assessed 
using the quiescent background event files produced by 
the ACIS team based on blank sky observations. These 
event files were used to simulate 20 200ks ACIS-I blank 
fields to match our data. The simulated images will give 
an estimate of the number of statistical spurious sources 
expected in 200ks data but will not account for pos- 
sible systematic spurious sources that may result from 
the more complex AEGIS-X data with its multiple over- 
lapping observations. Running our full source detection 
procedure on the simulated images we find there to be 
0.58 spurious sources detected per 200ks field per band, 
similar to the 0.5 sources per band estimated by N05. 
Therefore we expect there to be a total of no more than 
19 spurious sources in our catalog, equivalent to a con- 
tamination rate of less than 1.5%. 

4.6. Final AEGIS-X catalog 

The final AEGIS-X point source catalog is presented 
in Tables [5] and [TUl Table[5|gives the basic source proper- 
ties such as source position, source counts, and details of 
the detection significance. Column (1) give the source ID 
and column (2) gives the field ID, which identifies the in- 
dividual field in which the source was detected. Sources 
are listed in order of increasing RA. Column (3) gives 
the CXO object name. Columns (4) and (5) are the RA 
and Dec of the X-ray sources, respectively. These coordi- 
nates have been shifted to agree with the SDSS reference 
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Fig. 7. — (a) The residual offset in RA and Declination be- 
tween the matched Chandra and DEEP2 optical sources in the 
final AEGIS-X catalog, after the positional offsets were applied to 
each field (Table|3]l. (b) Total positional offset between the AEGIS- 
X and DEEP2 sources vs. off-axis angle. The median and RMS 
positional offsets are 0'.'40 and 0'.'79, respectively. 

frame of the DEEP2 images (Tabled. Column (6) is the 
positional error in arcsec, as detailed above. The source 
OAA, in arcmin, is given in column (7). Columns (8)~ 
(15) give the total source counts {N) and background 
counts (B) in the four analysis bands. Counts are given 
regardless of whether a source was detected in the band. 
Column (16) lists the bands in which the source is de- 
tected with Poisson probability < 4 x 10~^, where the 
bands are full (f), soft (s), hard (h) and ultrahard (u). 
Column (17) gives the lowest false detection probability 
found for the four bands; probabilities lower than 10~* 
are quoted in the table as 10~^. 

Table [TU] lists the flux and HR information for the 
sources. Column (1) again gives the source ID. Columns 
(2)-(5) give the fluxes calculated using the Bayesian 
method described in WS.'Sl which corrects for the Ed- 
dington bias. The errors are la values and the limits 
are Icr, or 68%, upper limits. In three cases (egs_0367, 
egs_0457, and egs_0633) sources that were significantly 
detected in a band had best estimate fluxes of zero, us- 
ing the Bayesian method. This discrepancy arises due to 



the different EEFs used for the photometry (90% EEF) 
and for the detection (70% EEF). In these cases we use 
source and background counts corresponding to the 70% 
EEF area to calculate the fluxes in the detected bands. 
Columns (6)-(9) give the observed frame fluxes in the 
four bands, where the fluxes and errors have been cal- 
culated using the classical method. The errors are la 
values, however here the limits are 99% upper limits. All 
the fluxes listed are in units of 10~^^ erg cm~^ s~^ and 
have not been corrected for Galactic or intrinsic source 
absorption. To correct for Galactic absorption the full 
and soft band fluxes should be increased by 2% and 4.2%, 
respectively. The corrections to the hard and ultrahard 
fluxes are negligable. Columns (10) and (11) give the 
Bayesian and classical HRs, respectively. The errors for 
the Bayesian HRs are la. For sources only detected in 
the full band, and with upper limits in both hard and soft 
bands, the clasical hardness ratio cannot be determined 
and is set to —99. Column (12) is a flag for the quality 
of the photometry. A flag of '1' indicates the presence 
of a nearby source that may be contaminating the pho- 
tometry. A flag of '2' indicates that another source was 
detected with the 90% EEF and that the photometry is 
likely heavily contaminated and the source position un- 
certain (j jS.Sp . All other sources have a flag of '0'. 

The AEGIS-X source catalog is publicly available in 
FITS table format at 



|http : / / astro . imperial . ac.uk/research/aegis/ cats . shtml 

The FITS version of the catalog includes more informa- 
tion about the sources, such as the Poisson probability, 
counts, PSF size, and exposure value in each of the 
bands. Postage stamp images of the sources are also 
available. 

5. OPTICAL AND IR IDENTIFICATIONS 

We identify secure optical and IR counterparts to 
the Chandra X-ray sources using the maximum likeli- 
hood (LR) method (|Ciliegi et al.ll2003l . l2005HBrusa et all 
l2007f ). This is the ratio between the probability that a 
source, at a given distance from the X-ray position and 
with a given optical or IR magnitude, is the true ID and 
the probability that this source is a spurious alignment. 
The LR ratio takes into account the X-ray, optical or 
IR positional uncertainty, the background density of op- 
tical (or IR) galaxies and the a priori probability of a 
counterpart with a given magnitude. 

For the positional accuracy of the X-ray sources we 
adopt a Gaussian distribution with standard deviation 
estimated as a function of off-axis angle and total num- 
ber of counts from the MARX simulations described in 
^ 34.41 The optical identification is a two pass procedure. 
In the first pass we estimate the a priori probability that 
an X-ray source has a counterpart of given magnitude in 
the input catalog (optical or IR). We use a fixed search 
radius of 2" to identify the X-ray sources with counter- 
parts in the input catalog. A total of 100 mock catalogs 
are then constructed by randomising the source positions 
of the input catalog and the cross-matching is repeated. 
The magnitude distribution of the counterparts in the 
mock and the real catalogs are subtracted to determine 
the magnitude distribution of the real associations. This 
is then used as a prior in the LR ratio estimation to ac- 
count for the distribution of X-ray sources in magnitudes 
and the total expected number of true associations. In 
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Fig. 8. — 0.5—2 keV flux versus ij-magnitude for sources with 
DEEP2 or CFHTLS counterparts. For X-ray sources with DEEP2 
counterparts we plot the DEEP2 magnitudes, for sources not 
detected in DEEP2 but detected in the CFHTLS then we use 
CFHTLS r' j^B converted to Rab (as described in i|5.1l l. oth- 
erwise we plot the DEEP2 upper limit. Lines of constant X- 
ray-to-optical flux ratio are plotted, according to the relation of 
[Hornschemeicr et al. (2001). The symbol shapes and colors are de- 
fined by the positions on the /x//opt plot as follows: (red squares) 
log(/x//opt) > 1, (black circles) -1 <log(/x//opt) < 1, (green 
inverted triangles) —2 <log(/x//opt) < —1, and (blue triangles) 
log(/x//opt) < -2. 

the second pass a search radius of 4" (equivalent to the 
3(7 positional uncertainty at large OAA) is adopted to 
identify all possible counterparts to X-ray sources and 
to estimate the LR of each one. Where more than one 
counterpart is found within the search radius we iden- 
tify the one with the highest LR to be the most likely 
counterpart. 

5.1. Optical identifications 

There is extensive optical photometric coverage of the 
AEGIS field from a number of telescopes, including 
SDSS, K PNO, Subaru and MMT (further details can be 
found in lDavis et al.l [20071 ). Here we consider two sets 
of optical photometric data, both taken at the Canada- 
France-Hawaii Telescope (CFHT), as our primary match- 
ing datasets. The DEEP2 phot ometric imaging i n the B, 
R, and / bands is described in lCoil et al.l (j2004D . These 
data almost the entire Chandra AEGIS field and are com- 
plete to a limiting magnitude of Rab = 24.1. The astro- 
metric accuracy of the photometric catalog is estimated 
to be G'.'2. 

We also consider deeper CFHT data that were ob- 
tained for a 1 deg^ area as part of the CFHT Legacy 
surveyi° (CFHTLS). The CFHTLS includes data from 
five filters: m*, g' , r' , i' , and z', and are complete to 
i'ab = 27.0. The overlap between the CFHTLS imaging 
area and the Chandra area is 0.33 deg^, considerably 

see |http://www.cf ht.hawaii.edu/S cience/CFHLS/| and 
|http: //www. ast.obs-mip.tr/article204. html 



smaller than the overlap with the DEEP2 optical imag- 
ing described above. The CFHTLS catalogs we are using 
are from the T0003 release, in which the astrometry was 
calculated using the USNO-Bl catalog. The astrometric 
uncertainty of the CFHTLS catalog is estimated to 
be about 0'.'3. Details about the data and catalog 
are giv en the TERAPIX pages relevant to that release: 
|http : //ter apix . iap . f r/article .php?id^rticle=556' 

Using the maximum likelihood method we search for 
optical counterparts to the AEGIS-X sources in the 
DEEP2 survey using the DEEP2 i?-selected catalog 
to estimate the surface density of background optical 
sources, R being the deepest band in the DEEP2 survey. 
The i' band is the deepest band in the CFHTLS survey 
so we use the i'-selected catalog to estimate the surface 
density of background optical sources for the CFHTLS. 
For both the DEEP2 and CFHTLS surveys we consider 
counterparts with LR>0.5 as secure, giving a '^6% con- 
tamination rate in both cases. 

Of the 1317 ^i?G/5'-X sources that are covered by the 
DEEP2 survey, which is complete to Rab = 24.1, we find 
that 897 have secure R band counterparts an optical ID 
rate of 68.1%. The deeper CFHTLS data, complete to 
i'AB — 27.0, overlaps with 703 AEGIS-X sources, and 
we find that 578 (82.2%) of them have secure i' band 
selected optical counterparts. Taking into account both 
catalogs, 1013 of the 1325 AECIS-X somces have optical 
counterparts, with 13.95 ^ Rab ^ 26.65. At the limit to 
which all of the optical coverage is complete, Rab — 24.1, 
we find an optical counterpart rate of 57%. The opti- 
cal ID rates found for the AEGIS-X sources are similar 
to those found for other dee p Chandra surveys. For in- 
stance, in the deeper CDF-N iBarger et aD (j2003f ) found 
optical counterparts for 85% of the sources to i? ^ 26.4, 
with an estimated 15% contamination rate. Considering 
only R ^ 24.0 sources, they found an optical ID rate of 
55%, in agreement with the ^i?G/5'-X findings. 

Table [TT] gives the results of the optical counterpart 
analysis. The table lists the R magnitude for every source 
with a secure DEEP2 counterpart, along with the offset 
between the X-ray and DEEP2 source (Sopt-x), the LR, 
and the DEEP2 counterpart ID number. The i' mag- 
nitude, offset between the X-ray and CFHTLS source, 
LR, and counterpart ID number are also given for every 
secure CFHLTS counterpart. 

Figure [8] shows soft X-ray flux versus R magnitude for 
the sample, for all the sources with secure optical coun- 
terparts. For X-ray sources with DEEP2 counterparts we 
use the DEEP2 R magnitude. For sources not detected in 
the DEEP2 survey, but securely detected in the CFHTLS 
then we convert the CFHTLS r' ab magnitudes to Rab 
using 

R = r' - 0.0576 - 0.3718[(/ - i') ~ 0.2589] (6) 

(jBlanton &: Roweii I2007T). This transformation is cor- 
rect for the iBessell et al] (I1998D R fiher which differs 
slightly from the CFHT fiher used in the DEEP2 sur- 
vey, but should be correct to within 0.1 magnitudes 
or less. Upper limits are plotted, where appropriate. 
The ratio of X-ray to optical flux is a useful diagnos- 
tic for the determining the nature of X-ray sources. 
Over several decades in flux classic AGN, both nar- 
row and broad line, exhibit X-ray to optical flux ratios 
of -1 <\og{fx/fopt) < +1 (e.g. iSchmidt eTldl II998I: 
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lAkivama efall 120001 : [Alexander et all l2001h . while at 
lower X-ray to optical flux ratios low luminosity AGN, 
starburst galaxies, a,nd normal galaxies emerge (e.g. 
[Alexander eraIll200H [larger et al.ll2002f) . The sources 
in the AEGIS-X catalog cover over four orders of mag- 
nitudes in fx-^o-fopt ratio. As expected, the majority of 
the AEGIS-X sources populate the AGN region of the di- 
agram but a significant number of sources fall within the 
regions normally populated with starburst or composite 
star-forming/ AGN galaxies (—2 <loK{.fx/ fop t) < ~1; 



e.g. 



-forming/ AGN galaxies (-2 <log(/x//op t) < 
[Bauer et al.ii2002t .Alexander et al.|[2002D and 



mal galaxies {log{fx/fopt) 
[2001 . 



< 



nor- 



2: [Hornschemeier et al.l 



5.2. IRAC identifications 

The majority of the AEGIS field is cover ed by sensitive 
Spitze r/IRA C data, described in d etail bv [Barmbv et al.1 
(|2006( 1 and [Barmbv et all (|2008f l. Briefly, these data 
cover the whole length of the strip but with a width 
of only 10 ^ Here we consider only the 3.6/im-selected 
catalog of iBarmbv et al.l (|2008f ). which has a 50% point 
source completeness limit of 0.9 fiJy at 3.6 fim (or 
'7i3.6/im(AB)=23.8). The astrometric uncertainty of the 
IRAC catalog is estimated to be 0'.'35. A total of 882 
Chandra sources are covered by the Spitzer region. Us- 
ing the maximum likelihood method, and again consider- 
ing only matches with LR>0.5, we find that 830 (94.1%) 
have secure 3.6/im counterparts. This cut on LR gives a 
^--^1% spurious counterpart rate. The 94% IRAC ID rate 
found here is consistent with that fo und in the 250ks 
Exten ded Chandra Deep Field-South; [Cardamone et all 
(|2008f) found IRAC counterparts for 90% of the X-ray 
sources to an IRAC flux limit of 0.63 ^Jy at 3.6 ^m. 

Table [11] lists the 3.6/im magnitude for each source 
with a secure IRAC counterpart, along with the off- 
set between the X-ray and IRAC source (Sm-x), the 
LR, and the IRAC ID from [Barmbv et al.[ ([2008). Fig- 
ure [9| shows soft X-ray flux versus 3.6/im magnitude for 
the sample, for the sources with an IRAC counterpart. 
Galaxies identified as classic AGN in Fig [5] (black circles 
and red squares) occupy the region of the plot defined by 
-2 <log(/x///fl) < 0. The sources in the AEGIS-X cat- 
alog cover over three orders of magnitudes in fx-to-fiR 
ratio. The AEGIS-X sources without an optical coun- 
terpart are on average fainter at 3.6/im and exhibit a 
smaller range in fx-to-fm ratio than those sources with 
an optical counterpart. 

6. CONCLUSIONS 

We have presented a catalog of X-ray point sources 
from the 0.67 deg^ AEGIS-X survey. Using a detec- 
tion threshold of Poisson probability < 4 x 10~^ we de- 
tected 1325 independent point sources down to on-axis 
flux limits of 2.37 x lO^^^ (0.5-10 keV), 5.31 x 10"" 
(0.5-2 keV) and 3.76 x lO'^^ erg cm'^ s^^ (2-10 keV). 
For each source we determine the X-ray flux in four 
bands (full, soft, hard, and ultrahard) using a Bayesian 
method that corrects for the Eddington bias, in addition 
to the standard flux calculation method. Optical and IR 
counterparts to the X-ray sources were identified from 
the DEEP2, CFHTLS, and Spitzer/lKkC surveys of the 
EGS, and the results presented. We find that 76% and 
94% of the AEGIS-X sources have secure optical and IR 
counterparts, respectively. 
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Fig. 9. — 0.5-2 keV flux versus IRAC 3.6/im magnitudes for 
sources with secure IRAC counterparts. Lines of constant X-ray- 
to-IR flux ratio are plotted, using the 3.6/im zero-points given by 
IReach et aT] ll2005l '). (Top) Sources with secure optical counter- 
parts: symbol shapes and colors are the same as Fig |S] (Bot- 
tom) Sources without secure optical counterparts. A few bright, 
'"S.e/tm < 19-0, sources that do not have secure optical counter- 
parts listed in Table 1111 have been removed from this plot as vi- 
sual inspection revealed clear bright, optical counterparts with mis- 
matching optical positions, such as saturated stars. Such sources 
are commented in Table [TT] 

A number of the analysis methods used in this work 
were specifically developed for Chandra surveys such as 
this one. The source catalog was produced using our 
own source detection procedure that is described in de- 
tail. Simulations verified the validity of this detection 
procedure and predict a very small, < 1.5%, contamina- 
tion rate from spurious sources. The photometry method 
employed corrects for the Eddington bias. A method for 
calculating accurate PSFs for multi-observation datasets, 
such as the AEGIS survey, is also presented. Sensi- 
tivity maps and sensit ivity curves were mad e following 
the new method of Ge orgakakis et aO ()200^. which ef- 
ficiently and correctly accounts for the observational bi- 
ases that affect the probability of detecting a source at a 
given flux. 



AEGIS-X 



15 



Finally, all of the data products described in this 
paper, including reduced event files, images, exposure 
maps, sensitivity maps, sensitivity curves, and PSF 
tables (which can be used for any Chandra ACIS-I data) , 
are available via the public website 
http : / / astro . imperial . ac . uk/r es earch/ aegis/[ 
The AEGIS-X catalog and the optical and IR counter- 
parts catalog are also made available. 



We thank the referee for helpful comments that im- 
proved clarity of the paper. We thank those who 
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TABLE 9 

Chandra AEGIS-X SOURCE CATALOG: BASIC SOURCE PROPERTIES 



05 









RA 


Dec 






FB cts 


SB cts 


HB cts 


UB cts 


Detection 




ID 


Field ID 


CXOEGS 


(J2000) 


(J2000) 


Pos. err 


OAA 


N 


B 


N 


B 


N 


B 


N 


B 


bands 


Pmin 


/I ^ 

(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


egs.OOOl 


egs7.198 


J141410. 9+520503 


213.545293 


52.084321 


1.33 


11.68 


28 


12.6 


14 


3.0 


15 


9.9 


8 


6.6 


s 


10-6.75 


egs_0002 


egs7.129 


J141418. 2+520856 


213.575718 


52.149002 


1.33 


10.11 


49 


21.0 


24 


5.4 


26 


17.2 


13 


10.6 


fs 


10-8.00 


egs_0003 


egs7.111 


J141419. 0+520404 


213.579288 


52.067991 


1.33 


10.82 


59 


12.0 


31 


2.9 


29 


9.2 


11 


5.6 


fsh 


10-8.00 


egs_0004 


egs7.128 


J141423. 0+520809 


213.595988 


52.135866 


1.33 


9.35 


49 


24.4 


21 


6.1 


32 


19.2 


16 


12.3 


fs 


10-6.32 


egs_0005 


egs7.193 


J141423. 2+520708 


213.596674 


52.118981 


1.33 


9.39 


38 


23.5 


16 


5.4 


22 


18.9 


14 


12.7 


s 


10-5.41 


egs_0006 


egs7.183 


J141424.0+520137 


213.599946 


52.027186 


1.33 


11.36 


42 


16.0 


18 


4.4 


23 


11.8 


13 


7.3 


f 


10-8.00 


egs_0007 


egs7.034 


J141426. 5+521136 


213.610293 


52.193375 


1.33 


9.41 


60 


18.7 


25 


4.4 


35 


15.4 


10 


10.1 


fsh 


10-8.00 


egs_0008 


egs7.110 


J141428. 2+520346 


213.617692 


52.062882 


1.33 


9.67 


58 


16.6 


39 


4.3 


19 


13.1 


5 


8.5 


fs 


10-8.00 


egs_0009 


egs7.122 


J141429. 8+520659 


213.624364 


52.116556 


0.87 


8.40 


72 


21.6 


37 


5.2 


35 


17.8 


15 


11.8 


fs 


10-8.00 


egs_0010 


egs7.127 


J141431. 1+520357 


213.629634 


52.065995 


1.33 


9.19 


53 


22.6 


31 


5.7 


25 


17.6 


12 


11.3 


fs 


10-8.00 


Note. — TablG[9]is publi 


.shed in its entirety in the 


electronic edition 


of tho journal. 



























TABLE 10 

Chandra AEGIS-X SOURCE CATALOG: SOURCE fluxes and HRs 

= H 

Bayesian flux Classical flux Bayesian Classical Phot. „ 



ID 


/0.5-10 




/0.5-2 


/2-IO 


h 


3-10 


/0.5-10 




/0.5-2 


/2-IO 




i-10 




HR 




HR 


flag 


(1) 


(2) 




(3) 


(4) 




(5) 


(6) 




(7) 


(8) 




(9) 




(10) 




(11) 


(12) 


egs_0001 


5.I6+I13 


1, 


f,o+0.73 
•"0-0. 60 


< 3.74 


< 


3.93 


"•4^-2.19 


2, 


nQ+0.89 

■'-'•-'-0.68 


< 6.70 


< 


9.52 


-0. 


.37t[;' 


,30 
,32 


-1.00 





egs.0002 


A.SOtlfr 


1 


qn+0.43 
•■-"'-0.37 


< 2.79 


< 


2.31 


c- 00 + 1.53 
1 33 


1, 


cc;+0.50 
■""^-0.41 


< 3.74 


< 


5.51 


-0, 


.381°: 


,28 
,24 


-1.00 





egs_0003 


17.68^3-11 


4, 


c:7+0.99 
•^'-0.88 




< 


5.79 


i8-42+i:^g 


4, 


OC- + 1.15 
■°"'-0.96 


13.99t377 


< 


8.90 


-0. 


.22l°' 


,18 
,14 


-0.22 





egs_0004 


r. 70 + 0.95 

^- '■-'-0.87 





71+0.27 
•'1-0.23 


1 70+1. 08 
1^ '^-1.72 


< 


1.82 


icr + 1.03 
•-'•l>'-0.89 


0, 


00+0.32 
■°'-'-0.25 


2-99tl:|l 


< 


3.81 


-0, 


,131° 


,28 
,25 


-0.12 





egs_0005 


1 10+0.57 
1-1^-0.82 


0, 


44+0.23 
•^^-0.20 


< 1.15 


< 


1.49 


1 70+0.90 
!• '^-0.76 


0, 


^•7+0.28 
■"'-0.21 


< 2.73 


< 


3.92 


-0. 


.4ois;' 


,30 
,40 


-1.00 





egs.0006 


9.68t2;80 


2 


in+0.82 
•1^-0.69 


5 04+2. 80 
■^•■^^-3.70 


< 


6.31 


10.72t^;ii 


2, 


40+0.97 


00+4.31 

^•^■J-3.49 


< 


9.69 


-0, 


■i6lS: 


,27 
,25 


-0.15 





egs.0007 




1, 


f,o+0.45 




< 


1.95 




1, 


7f.+0.52 
■ '"-0.42 


07+2.49 
°-'^'-2.10 


< 


5.33 


-0. 


.07l°' 


,20 
,18 


-0.07 





egs_0008 


7 qq+l-SS 

i.yo_i.46 


2 


Q9+0.57 
•^^-0.51 


< 2.12 


< 


1.59 


8.36+1-53 


3 


07+0.65 
■^^'-0.55 


< 3.62 


< 


5.28 


-0, 


,691° 


,14 
,18 


-1.00 





egs_0009 


4.00_Q 7g 


1, 


00+0.26 
•^■^-0.23 


2 ^9+11" 
^•■5^-1.02 


< 


1.27 


4 7fi+o:9o 

^- '"-0.80 


1, 


qn+0.29 
■■^"-0.25 


^•='='-1.03 


< 


2.67 


-0. 


.331° 


,18 
,17 


-0.33 





egs_0010 


1C+O.86 
■J-l'J-0.79 


1 


it:+0.29 
•l''-0.25 


< 1.47 


< 


1.26 


4(5+0.95 
•^•™-0.82 


1, 


ncr+0.33 
■^"-0.27 


< 2.39 


< 


3.42 


-0, 


,541°: 


,21 
,23 


-1.00 






Note. — Tabic flO] is published in its entirety in the electronic edition of the journal. 



TABLE 11 

Optical and IR counterparts to the AEGIS-X sources 







DEEP2 




CFHTLS 






IRAC 




ID 






LRopt 


ID 


m.i^ Sopt-x° LRopt ID 




R — X ^ 


LR/fl 


ID" 




("AB) 


(") 






CAB1 (") 


(AB^ 


{"^ 






egs.OOOl 


22.84 


2.02 


1.55 


11021738 












PCTQ 0009 


zo.oo 


1.04 


4.15 


11028286 












egs_0003 


21.11 


1.65 


7.14 


11015296 












egs_0004 


22.00 


1.49 


7.68 


11028607 












egs_0005 


21.80 


1.15 


13.01 


11021753 












egs_0006 


>24.1 










>23.8 








egs_0007 


22.30 


0.64 


31.25 


11035367 












egs_0008 


>24.1 










20.82 


1.31 


9.46 


J141428.13+520347.1 


egs_0009 


23.94 


0.70 


10.04 


11021149 












egs.OOlO 


22.66 


1.26 


6.69 


11014844 




20.25 


1.09 


25.36 


J141431. 06+520358. 6 


Note. — 


Tabic fm is 


published 


in its entirety in the electronic edition of the journal. 











^ If optical or IR coverage exists but no secure counterpart was detected then an upper limit for the magnitude is given. X-ray sources with no optical or 
IR coverage are denoted as '. . . '.'^ The offset between the X-ray and optical position.*^ The offset between the X-ray and IRAC position.'^ Truncated IRAC 
object ID number from IBarmbv et aT] l l200gD . Object IDs should be preceded by EGSIRAC. 



